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Application of the GRC Stirling Convertor 
System Dynamic Model 


Timothy F. Regan and Edward J. Lewandowski 
Sest, Inc. 

Middleburg Heights, Ohio 44 1 30 


The GRC Stilling Convertor System Dynamic Model (SDM) has been developed to 
simulate dynamic performance of power systems incorporating free-piston Stirling 
convertors. This paper discusses its use in evaluating system dynamics and other systems 
concerns. Detailed examples are provided showing the use of the model in evaluation of off- 
nominal operating conditions. The many degrees of freedom in both the mechanical and 
electrical domains inherent in the Stirling convertor and the nonlinear dynamics make 
simulation an attractive analysis tool in conjunction with classical analysis. Application of 
SDM in studying the relationship of the size of the resonant circuit quality factor (commonly 
referred to as Q) in the various resonant mechanical and electrical sub-systems is discussed. 


Nomenclature 

AreaD Area of displacer 

m p , nid moving mass of convertor piston (p) or displacer (d) 

Qp> Qd quality factor of piston’s (p) or displacer’s (d) resonant circuit 
s pole location in the complex plane 

SdP sum of pressure drops across heater, regenerator and cooler 

X p , X d amplitude of piston (p) or displacer (d) 

co 0 operating frequency 


I. Introduction 

The free-piston Stirling convertor is being developed as an efficient and reliable power source for NASA 
exploration missions. The 110 watt Stirling Radioisotope Generator (SRG110) project is being conducted by the 
U.S. Department of Energy with system integrator Lockheed Martin, Valley Forge, Pennsylvania. The project will 
integrate two Stirling convertor assemblies (SCA’s) into a space power generation system. 1 The SCA’s used are the 
55 We Technology Demonstration Convertor (TDC) 2 manufactured by Stirling Technology Company (STC), 
Kennewick, Washington. Predictions about how the SRG110 will behave under mission conditions can be made by 
test and simulation. The simulation tool under development at NASA’s Glenn Research Center (GRC) is the Stirling 
Convertor System Dynamic Model (SDM). It can be used to predict SRG110 dynamic performance under mission 
conditions of interest. The simulation results can be reviewed and decisions can be made about follow-on studies or 
hardware tests based on the simulations. The model has been configured to take thermal power as input and give 
electrical power, temperatures and mechanical motion as outputs. In this paper, use of the model is demonstrated by 
presenting simulation results and comparing them to test data. Development of the model is treated in Ref. 3 through 
5. Another application of the model is as an aid in controller design. This subject was treated in Ref. 6 and 7 and is 
not discussed further here. The final application that is discussed is the use of simulation as an alternative to 
classical analysis techniques in understanding the behavior of a convertor systems comprised of more than one SCA. 
Berchowitz has presented full details concerning the dynamics of a single SCA using classical analysis. 8 Simulation 
can be a more practical analysis technique than classical analysis to answer questions about electronic controllers 
and multiple SCA interactions. 
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II. Application: Prediction of Off-nominal Behavior 

A number of simulations have been performed to date as part of the SDM validation effort. Reviewing these 
simulations gives some insight into the model’s application as a means to investigate system behavior under off- 
nominal conditions. This is a step in evaluating the system reliability. In the reliability evaluation, component values 
and operating conditions are made to vary within a band of values. As they vary, it is to be expected that at some 
point, operation will be affected. For example, if the cold-end temperature were increased by a certain percentage, 
one would expect the convertor performance to be reduced. A computer code like Sage 9 can be used to evaluate 
losses and output power at the off-nominal cold-end temperature — it was written to evaluate Stirling machine 
performance. The SDM is able to simulate system dynamics under the same off-nominal conditions. What can look 
like reduced but acceptable performance under a Sage 9 analysis can conceivably register as a failure in a reliability 
study when dynamic performance analyzed by SDM is considered. The reason is that the imbalances which cause 
system instability can be detected by SDM. By analyzing unbalanced operation cases with both Sage 9 and SDM, the 
boundaries between acceptable performance and unacceptable performance can be mapped not only for cold-end 
temperature variation but for any parameter or combination of parameters desired. Table 1 shows some descriptions 
of simulations that have been run. 


Table 1. — Simulations performed as part of SDM validation. 


Simulation 

Simulation Description 

System 

A 

Single simulation, changing values of SCA-1 Tuning capacitance. 

Single SCA 

B 

Multisimulation, SCA-1 Tuning capacitance as parameter. 

Dual-opposed 

C 

Multisimulation, SCA-1 and -2 Tuning capacitance as parameter. 

Dual-opposed 

D 

Multisimulation, SCA-1 cold end temperature as a parameter. 

Dual-opposed 



Figure la. — Dual opposed Stirling convertors 
flexibly mounted to a rigid base. 






SCA1 


Controller 


\\\\\\ 


Figure lb. — Single convertor rigidly 
mounted to a base. 


The single system listed in Table 1 includes one SCA with its case mounted to ground. The single system is 
shown in Figure lb. Simulations listed in Table 1 as dual systems include two SCA’s in the general configuration 
shown in Figure la. The mounting mass, stiffness, and damping shown in Figure la were chosen to match the 
mounting conditions measured on the test stands of the GRC Stirling Research Lab. In the space power system 
evaluation, the mounting must reflect the vehicle’s mounting conditions. In validation, the exact test conditions are 
modeled because the results will be compared to test results. 
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A. Simulation of tuning capacitor changes 

Simulation A is a single SCA simulation using off-nominal values of tuning capacitance. The tuning capacitor is 
a system component used to correct the power factor of the alternator so that the alternator current can be 
minimized. There are mass reduction benefits associated with minimizing alternator current. The tuning capacitor 
forms a resonant circuit with the alternator stator inductance. The tuning capacitance required to meet the minimum 
stator current design objective also makes the resonance frequency of the electric circuit the same as the operating 
frequency of the convertor. The value of the tuning capacitor therefore can have an effect on system dynamics. 

In the simulation, a fixed input power was applied to the hot end. The hot-end temperature was allowed to float. 





Output Power 
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Figure 2. — Comparison of test result with simulation result — output power (left) operating frequency (center) 
and terminal voltage (right) Configuration of Figure lb was used with off-nominal tuning capacitors. 

Tuning capacitance shown as abscissa. 

The cold-end temperature was fixed. The capacitance value was increased above and below the design value to 
match the test cases that were run. The simulated power, frequency and terminal voltage tracked the values 
measured in test in general. In testing, tuning capacitance was increased to a factor of 2.16 times nominal. Over this 
range, no instability was observed. The graphs of Figure 2 track the performance of both the test and the simulation 
in relation to the value of tuning capacitance connected. The alternator terminal voltage and operating frequency 
were predicted well by the model. The output power calculated by SDM is high because of the isothermal 
assumption used to model the thermodynamic cycle. SDM does not calculate all the known losses because the 
development plan calls for SDM to run Sage in background to improve the thermodynamic modeling of the cycle. 
Sage is a program that calculates Stirling machine performance including losses. When the Sage interface upgrade is 
accomplished, the output power simulated will be more accurate. The simulated displacer phase angle although not 
shown in Figure 2 was also low (by about 11°.) The phase angle discrepancy is also an artifact of the isothermal 
assumption and will be also corrected with the Sage interface. 

Values of tuning capacitance larger than what is shown in Figure 2 were also used in the single SCA simulation 
in order to predict the value at which instability will occur. A tuning capacitance value of 50 pF caused instability in 
the simulation. In that case, the cause of instability was that the load limit of the controller had been reached. It is the 
controller that determines the alternator current that in turn provides the electromagnetic damping force which 
extracts real power from the piston. The controller accomplishes this by reducing or increasing load resistance seen 
by the alternator until the circuit voltage satisfies a set point. With the proper value of tuning capacitance connected, 
the stator current required to extract rated power from the piston is well within the capability of the controller. In the 
case of an off-nominal value of tuning capacitance the power factor is lower and more current is required. With a 
50 pF tuning capacitor, the controller could not synthesize the electrical resistance required for the necessary 
current. The controller can synthesize any load required by normal operation of the convertor plus a margin. At 
some point in off-nominal operation, the margins will be exceeded and in the case of the single SCA simulation, the 
margin was exceeded at 50 pF of tuning capacitance. The value has not been tested in hardware and in fact it is 
likely to experience instability with a tuning capacitor lower than 50 pF. 

As has been pointed out, one cause of instability of a single SCA when operated with off-nominal tuning 
capacitors is the limit of load available from the controller. Another is a limitation of flux flowing in the alternator 
core. The alternator is modeled in the SDM with infinite flux-producing capability and so the present version of the 
SDM does not model the instability caused by flux limitation. The importance of mass reduction in aerospace 
applications implies that alternator cores, which are made of iron, must be sized to carry design flux plus margin and 
no more. For that reason, it is important for the SDM to incorporate flux limits in its alternator model. 
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The modeling platform used by SDM has a magnetics physical domain available. The magnetics physical domain 
contains all the circuit elements necessary to model the flux limitation of the core used in an SCA. This suggests 
another upgrade to the SDM: incorporation of a non-linear alternator flux versus excitation characteristic. 

B. Simulation of unbalanced tuning capacitance 

Simulation B is a simulation of a dual-opposed system at about half rated power in which the value of one of the 
tuning capacitors was changed while the other remained fixed. The tuning capacitance for SCA-1 was changed. That 
of SCA-2 remained unchanged. The heat input was fixed and the hot-end temperatures of each SCA was allowed to 




Figure 3. — Comparison of test with simulation results for SCA-1 and SCA-2 output power in dual-opposed 

simulation (Fig. lb) with unbalanced tuning capacitors. 




Figure 4. — Comparison of test with simulation results for SCA-1 and SCA-2 piston amplitude in dual- 
opposed simulation (Fig. lb) with unbalanced tuning capacitors. 


float. In this simulation, the unbalance was observed to cause motion of the case mass. The motion of the case mass 
in turn causes additional effects in performance. Simulation B was performed at power levels that were well below 
rated power. The flux limitation discussed under simulation A is related to the alternator current level. In simulation 
B, the alternator current level never exceeded the rated alternator current, due to the reduced power level at which 
the simulation was conducted. As a result, the flux limit of the hardware was not exceeded in the simulation. The 
simulation results are compared with the test results in Figures 3 and 4. Again, due to the isothermal assumption of 
SDM, the output powers in the simulation were greater than those of the test results but the trend with increasing 
tuning capacitance was in agreement with test results. The model predicts the dynamic performance as represented 
in Figure 3 by piston amplitudes. The frequency of operation, not shown tracked the test results 1 Hz higher at each 
of the 4 test points. 
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Figure 5. — Instability as a result of balanced tuning 
capacitors whose value exceeded the nominal value 
by factor greater than the margin in simulation C. 



Figure 6. — Terminal voltage of SCA1 in a dual 
simulation according to Figure la with 20uF 
tuning capacitors and 40w output. 13Hz 
subsynchronous oscillations shown. 


C. Simulation of balanced but off-nominal tuning 
capacitance 

Simulation C was performed to further investigate 
tuning capacitance variation. The dual-opposed 
arrangement of Figure la was again used. The tuning 
capacitors that are connected to each SCA were kept 
equal to each other but a value of capacitance 1.23 
times the design value was used. This was a value that 
was known to cause instability in test. The test 
demonstrated unstable operation of the system when 
output power was 40 We. The simulation used the 
constant input power level required to convert 40 We 
per SCA. The simulation predicted instability as shown 
by the exponential rise in case motion shown in 
Figure 5. 

Another prediction of the SDM was the form and 
frequency of the sub-synchronous oscillations observed 
in the testing. In test, the onset of instability was 
accompanied by a 9 Hz sub-synchronous beat 
frequency. One second’s worth of the alternator 
terminal voltage is shown in Figure 6. In simulation, a 
beat frequency of 13 Hz was predicted as shown in 
Figure 6. 

D. Simulation of Unbalanced Cold-End Temperature 

Simulation D simulated the effects of an unbalanced 
cold-end temperature on operation. The test and simulation 
are of practical interest because the cold end temperatures of 
an SRG110 may vary if one SCA is exposed to more 
sunlight than the other SCA. In such a case, the cold-end 
temperature and the bounce-space temperature may be 
presumed to be different. Three cases were considered in 
test and simulation. Simulation D considered balanced 
temperatures at the extremes and an unbalanced 
temperature in the center. In the unbalanced case, both the 
cold-end temperature and the bounce-space temperature are 
different. The test results are compared to simulation results 
in Figure 7. 



SCA2 Output Power 



Figure 7. — Comparison of test and simulation output power for dual SCA configuration according to 
Figure la with constant heat input and shade-to-sun temperature transition. 

Case 1 - full shade, Case 2 - transition, Case 3 - full sun. 
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SCA1 results show excellent agreement between test and simulation. The output power of SCA2 behaves 
differently in simulation than in test in the middle temperature case. The transition case is the case for which the 
SDM calculated most unbalance in the dynamics. 

III. Application: SDM as an Aid to Analysis of Convertor Systems 

In Ref. 8, Berchowitz has described an elegant treatment of the linear dynamics of free-piston machines. It used 
complex analysis to derive expressions for complete dynamics including starting conditions and operating frequency 
of a single convertor. The analysis used the calculated Q (or quality factor) of the piston and displacer resonant 
circuits to derive the pole locations in the complex plane as the operating point of the machine is varied. The linear 
analysis derives Q from linearized equations of piston and displacer motion. SDM can compute the value of Q from 
measurements of piston and displacer motion. The Q of a resonant circuit is the ratio of stored energy to dissipated 
energy per cycle. The mass, spring, and damper of the piston in a free-piston machine is such a circuit. Its kinetic 
energy can be calculated as a function of its frequency and amplitude as in Eq. (1), assuming sinusoidal motion. 

KE = ^--m P - co o 2 X p ( 1 ) 


The energy converted to electrical power is the energy dissipated per second, Pe. The energy dissipated per cycle 
is PJfo where f, is the operating frequency in Hz. Then the piston Q is calculated as the ratio of stored energy to 
dissipated energy per cycle. 


Qp 


m p ■ (0 2 X 2 

IPjfo 


( 2 ) 


The value of Q P is seen to change as the piston amplitude and output power change. SDM has been configured to 
calculate the value of Q P with time. Typical results are shown in Figure 4. 

A similar definition holds for the displacer mechanical circuit. The power dissipated in the displacer circuit is 
calculated by SDM as 


P D = SdP ■ AreaD ■ co o X D / 


( 3 ) 


where SdP is the sum of the pressure drops across the heat exchangers and regenerator and X D is the displacer 
amplitude. The displacer Q is then 


Qd 


m D ■ 0 ) o 2 X d 2 

2 P D /f„ 


( 4 ) 



Figure 8. — Location of the piston and displacer poles for 
several operating points as calculated by SDM. Arrows 
show the location as output power increases. 
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the piston and displacer Q values were calculated by SDM and used to determine pole locations in the complex 
plane. The calculation is from Eqs. 5 and 6 obtained directly from Ref. 8. The results are plotted in Figure 8. 


S dl ’ S d2 


= CO, 


-\l{A-7TQ d )±j. 



( 5 ) 


Spi,S p2 =CO p - 1 


-H(A-7T-Q )±j \l- 


4 -x-Q d 


(6) 


The piston poles are represented by the progression of squares moving to the left in Figure 8. The poles move to 
the left as output power increases. The displacer poles of the convertor analyzed move to the right as power 
increases. As pointed out by Berchowitz, the operating frequency of the final operating point of the convertor can be 
determined graphically from the piston and displacer pole locations. s Knowledge of the location of the piston and 
displacer poles in the complex plane can also be used to advantage in system identification. 

In addition, modeling of a system with SDM makes all pressures, pressure drops, motions and temperatures 
available for display and further analysis. Phasor diagrams are very useful tools to be used in visualizing operation 
of a Stirling convertor system. The SDM makes magnitudes and angles of all periodic signals available to the 
analyst so that phasor diagrams can be constructed. 


IV. Conclusion 

It has been shown that system behavior of interest in dynamics and reliability studies can be obtained from 
NASA GRC’s Stirling Convertor System Dynamic Model (SDM). The dynamic effects of unbalanced operation can 
not be evaluated using other software tools available. Examples considered included unbalanced temperatures and 
unbalanced tuning capacitance. Since the SDM is physics-based, it is also convenient to study unbalances in 
component characteristics - the regenerator fiber diameter as an example. 

Two important SDM improvement areas have been pointed out and the plan for incorporating them has been 
mentioned: 

1. ) Incorporation of more realistic modeling of the alternator core 

2. ) Incorporation of more sophisticated thermodynamic modeling of the cycle by running Sage in the background 

Use of a simulation model of a nonlinear aids understanding of system behavior by making all the states 
available for inspection and analysis. In a practical space application, it is necessary to understand the limits of 
parameter variation that cause unacceptable operation as in the reliability evaluation discussed. It is also necessary to 
understand which component values can be adjusted to obtain more reliable dynamic performance. 
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